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Abstract 

The behavior of a single impurity in a one-dimensional Luttinger liquid is nu- 
merically investigated by means of the density matrix renormalization group. 
By analyzing the finite size scaling behavior of the low energy spectrum, we 
confirm the theoretical prediction of Kane and Fisher [Phys. Rev. Lett. 68, 
1220 (1992)] both for attractive and repulsive interactions. Moreover, we cal- 
culate the exponent of the orthogonality catastrophe, which gives a further 

support to the above theoretical prediction. 
71.10.Pm,72.15Nj 
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The problem of an impurity in a one-dimensional (ID) interacting Fermi system (Lut- 
tinger liquid) has become quite popular in recent years for its implications to a variety of 
physical problems as for instance the behavior of quantum wires (see e.g. Refs. |1|,0]) or the 
tunneling through a constriction in the fractional quantum Hall regime 0]. 

The interesting feature of this problem lies in the fact that the competition between the 
impurity potential and the electron-electron interaction leads to quite surprising effects. If 
the interaction is repulsive, the electrons at low energy see the potential barrier as if it were 
effectively infinite. On the contrary, for an attractive interaction the effective scattering close 
to the Fermi energy is vanishingly small (see e.g. Ref. 0]). The two different behaviors have 
a quite simple physical interpretation: in ID the electron gas is quasi-ordered, i.e. it shows 
power law decaying correlation functions of the order parameter. Specifically, for repulsive 
interaction it has a charge density wave, while for attractive interaction a superconductive 
quasi-long-range order. It is then clear that the two different orderings lead to opposite 
behavior in the presence of a local potential, which tends to pin the charge density wave but 
has minor effects for a superconductor. 

/^From the theoretical point of view, the problem has many similarities with the Kondo 
effect. In fact the problem of a magnetic impurity in a metal can also be transformed into 
a one-dimensional model of electrons in the presence of a local potential. Moreover, in 
both cases a perturbation expansion in powers of the impurity potential breaks down at low 
energy due to the appearance of logarithmic singularities. A standard approach based on 
renormalization group leads in both cases to the conclusion that the impurity potential flows 
either to infinity or to zero, depending on the the sign of the exchange in the Kondo model 
or of the interaction in the Luttinger liquid. Therefore, analogously to the Kondo effect, a 
low energy description of an impurity in a Luttinger liquid needs a correct identification of 
the low energy fixed point. In particular, Kane and Fisher [^] (KF) argued that the fixed 
point corresponds to a chain disconnected at the impurity site for repulsive interaction. 

This interpretation has been recently questioned by Oreg and Finkel'stein |^ (OF), 
specifically in the context of the X-ray edge singularity. These authors find in fact an expo- 
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nent of the absorption spectra at the threshold which is not equal to the value compatible 
with a disconnected chain fixed point It is however unclear whether they would extend 
their criticism to all the results which have been obtained on the basis of the disconnected 
chain fixed point, or if they only doubt about the X-ray edge singularity results. In fact, for 
other properties, like for instance the finite size spectrum [^] or the conductance through the 
barrier there exist exact numerical calculations which support KF's hypothesis, while 
for the X-ray edge singularity only approximate analytical results are till now available P]. 

In this paper we investigate this problem by numerical methods, closely following the 
kind of analysis in Ref. [^. In particular we study the low energy spectrum for ID periodic 
chains of spinless fermions interacting either with attractive or repulsive interaction in the 
presence of a local potential. Our numerical results show that this spectrum is compatible 
with a disconnected chain fixed point for repulsive interaction and with a periodic chain if 
the interaction is attractive. We also calculate the exponent of the orthogonality catastrophe 
from the overlap between the ground states in the presence and in the absence of the local 
potential. This exponent is related to the X-ray edge singularity exponent, and therefore 
our results can solve the above mentioned controversy. We find that the orthogonality 
catastrophe exponent again points in favor of the KF hypothesis. 

The development of the density matrix renormalization group (DMRG) method has 
made possible the study of the ground state and low excited states of one-dimensional 
systems with large number of sites. Previous DMRG studies on 5 = 1/2 chains has proved 
also quite successful in determining the finite size behavior of the low excitation energies. 



correlations and Friedel's oscillation due to boundary effects |T0|. In this paper, we use 
DMRG to study the low energy spectrum for anisotropic S = 1/2 Heisenberg chains with 
a local perturbation. This model can be mapped through a Jordan- Wigner transformation 
onto a model of interacting spinless fermions with a local impurity potential. 

An S = 1/2 spin-anisotropic Heisenberg chain with open boundary conditions (OBC) is 
described by the Hamiltonian 
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where 5*^^ and S^^ are spin raising and lowering operators at site z, Sf being its z-component. 
For periodic boundary conditions (PBC) the L*^ and l'^* site are also coupled by the same 
exchange couplings J^y and J^. An impurity site potential is introduced by adding a nonzero 
local magnetic field at the 1** site in the periodic chain. An impurity bond is instead 
introduced by coupling the L*^ and 1^* site with exchanges J'^y and J'^ different from the 
bulk values Jxy and J^. 

Eggert and Affleck |^ have shown that the systems ffow to the OBC fixed point whenever 
the impurity bond strengths between the 1'^* and L^^ sites differs from the bulk values. They 
have studied odd length chains under global SU(2) symmetry {J^y = Jz and J'xy = J'z) by 
comparing the parity, multiplicity, and magnitude of energies of the ground state and the 
lowest excited states. We have extended their analysis by considering both even and odd 
length chains, either for site impurity or bond impurity, and also in the case when SU(2) 
symmetry is explicitly broken. 

Particularly, in the presence of spin anisotropy, we find numerically that the low energy 
spectrum in the presence of a local potential flows to the OBC spectrum if < Jz/Jxy < 1, 
and to the PBC spectrum if — 1 < Jz/Jxy < 0, thus in agreement with the theoretical 
prediction of Ref. [Q. For the sake of clarity, we are going to present our numerical data for 
the specific case of even length chains with bulk anisotropies Jz/Jxy = 0.5 and —0.5, which 
correspond, respectively, to repulsive and attractive interactions in the equivalent spinless 
fermion model. 

We introduce the impurity potential by bond couplings J'^ = bJ^ and j'^y = bJxy, and 
denote the anisotropy by a = Jz/ Jxy We calculate the ground state energy Eq and the first 
excited energy Ei for even length chains by DMRG. The ground state and the first excited 
state are the lowest energy states in the sectors with 2;-component of the total spin 5**°* = 
and 1, respectively. We use the parity of the ground state as a reference, and we consider as 
quantum numbers of the lowest energy states their 2;-component of the total spin and their 



parity. Each state is therefore represented by a fixed value of S^^~\ where + is the parity of 
the ground state and — denotes the opposite parity. The ground state is then represented 
by 0^ and the first excited state Ei by 1^. The two energy levels of the low energy spectrum 
we choose to study are 0~ and I"'', and we denote their energies by E2 and E^, respectively. 

For different impurity bond strength b, L[Ei{a,b) — -E'o(a, 6)] = irvp is a bulk property 
and therefore remains the same for fixed anisotropy a. This is shown in Fig.l for a = ±0.5 
and b = 0,0.1 and 1, (6 = and 1 are the values for OBC and PBC, respectively). The 
behavior of the low energy spectrum in the presence of the local potential is instead made 
evident by the fiow of the two energy levels E2 and E^. We use the energy unit nvp/L, i.e. 
we scale the energy as 



For attractive electron-electron interaction (a = —0.5), we plot in Fig. 2 e2(— 0.5,6) and 
e3(— 0.5,6) for different 6-values: 6 = 0, corresponding to OBC, 6=1 for PBC, and 6 = 0.1 
for an intermediate case. It is clear from Fig. 2 that e2(— 0.5,0.1) and e3(— 0.5,0.1) fiow to 
e2(— 0.5, 1) and e3(— 0.5, 1), respectively, i.e. they fiow to PBC. We just mention that other 
energy levels of the spectrum for attractive interaction — 1 < a < also fiow to those with 
PBC. For repulsive electron-electron interaction (a = 0.5), we plot in Fig. 3 62(0.5,6) and 
63(0.5,6) for different 6-values. The energy levels of the spectrum for repulsive interaction 
flow to those with OBC. 

The same conclusions can be drawn for other values of the spin-anisotropy a and other 
impurity bond strength 6, or by substituting the bond potential with a site potential provided 
by a local magnetic fleld, or else for the odd length chain case (which however has already 
been extensively discussed in Ref. [§). 

Therefore our numerical results show that PBC is the flxed point towards which the 
system in the presence of a local potential flows when the interaction is attractive, while 
OBC is the flxed point for systems with repulsive interaction. 



63(0,6) 



62(0,6) 



E2{a, 6) - Eoi^a, 6) 

Ei{a, b) — Eoia, b) ' 
Ei{a,b) - Eola,b) 

Ei{a, 6) — Eo{a, 6) 



(2) 
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An interesting quantity which can also be evaluated numerically and which further sup- 
ports the interpretation of the strong coupling fixed point in terms of a perfectly reflecting 



barrier is the exponent a of the orthogonality catastrophe |TT]]. This exponent is defined 
through the overlap between the ground state wave functions for a system of size L in the 
presence, and in the absence, |0o)) of the impurity potential. In particular. 



m^)-\-i) ■ (3) 

The orthogonality exponent is intimately related to the exponent of the X-ray edge singu- 
larity. Let us assume that the scattering potential is present if some localized level is empty 
and is absent otherwise. The localized state Green function can be shown to decay at long 
times like: 

Git) = {<Po\d^{t)d{OMo) 



where d is the operator which empties the localized level, Hq and H the Hamiltonians in 
the absence and in the presence of the scattering potential, respectively. In the absence of 
electron-electron interaction, the exponent a can be exactly determined |jlT|] and is given by: 



The phase shifts of the even {6e) and odd {6o) scattering channels can in turn be related to 
the properties of the barrier, since 

Se{o) = ^(^S+± tan"^ , (6) 

where the transmission amplitude t = \t\e^^+, and r is the reflection coefficient. 

In order to simplify the interpretation of the numerical results, we consider an anisotropic 
Heisenberg chain of length L. The ground state is in the subspace with total spin z- 
component Sz = 0. The impurity potential is provided by modifying the exchange of a 
single bond in the chain. This model maps onto a chain of interacting spinless fermions at 
half filling. The impurity potential is such as to preserve the particle-hole symmetry of the 
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Hamiltonian. As a consequence, the forward scattering phase shift 5+ = 0, since a finite 
5+ would imply a breaking of the particle-hole symmetry. Before presenting the numerical 
results, it is worthwhile to discuss what the theory would predict for the exponent a. As we 
already mentioned, according to Kane and Fisher the strong coupling fixed point towards 
which the model flows in the long wavelength limit consists of an Heisenberg open chain: the 
impurity has simply changed the boundary conditions from periodic to open. In this case 
the exponent is predicted to be a = 1/16 f^, which corresponds simply to take 5+ = and 
|r/t| — > oo in Eq.(P). As we previously said, the KF interpretation has been questioned by 
OF in the context of the X-ray edge singularity 0. Following them, we expand the Green 
function in powers of the impurity potential. Each term in the perturbation expansion 
can be evaluated exactly in the long-time limit and finally we find, in accordance with OF, 
that the imaginary time Green function G'(r) coincides with the partition function of a 
two-dimensional classical Coulomb gas confined on a line of length r. The relevance of the 
impurity potential translates in the language of the Coulomb gas into an increasing fugacity: 
the gas is therefore in the plasma phase. If r is much bigger than the screening length of 
the plasma Tgcr-, then OF have predicted that 

lim G(r) = e~^^''<>''^ const., (7) 

which would imply an exponent a = 0. As we see the two theoretical approaches lead to two 
completely different results; therefore the exact numerical approach we are going to describe 
turns out to be quite decisive. 

We denote the overlap integrals to be evaluated numerically and their exponents as: 

Op(L) = (0PBC|0) oc (1)", 

Oo{L) = {<PoBc\<P) oc il'f, 
where the ground states are \4>pbc) for PBC, \4>obc) for OBC, and |0) for chains with PBC 

but in the presence of a modified bond < 6 < 1. The exponents a and (3 are given by 

, lnOp(L + 2) -In Op(L) 

ala, o) = — — -7 , 

^ ^ lnL-ln(L + 2) .qs 

lnOo(L + 2)- In Oo(L) 
lnL-ln(L + 2) " 

7 



For anisotropy a = ±0.5 and impurity bond strength b = 0.1, 0.9, we calculate the ground 
states by DMRG method, and then calculate the exponents a{a,b) and (3{a,h). We plot in 
Fig. 4 the exponents as a function of — 2/[lnL + ln(L + 2)] for different sizes L. The figure 
shows that a{a, b) flows to 1/16 and /5(a, b) to for repulsive interaction (a = 0.5), and vice 
versa for attractive interaction (a = —0.5). 

The values a = and /5 = 1/16 which we find for attractive interaction are in agreement 
with the flowing of the system to the PBC fixed point. On the other hand, a = 1/16 and 
/3 = for repulsive interaction show that the system flows indeed to the OBC fixed point, 
thus in agreement with Kane and Fisher's interpretation of the strong coupling fixed point 
and with the analytical results of Refs. ||^ . 

In conclusion, we have investigated by numerical DMRG method the impurity problem in 
a Luttinger liquid. All our results, which include a detailed analysis of the finite size scaling 
behavior of the low energy spectrum and the evaluation of the orthogonality catastrophe 
exponent, suggest that at low energy the barrier acts as it were perfectly reflecting. The 
orthogonality catastrophe exponent result may help to clarify the controversy recently raised 
by the authors of Ref. ^ about the X-ray edge singularity |T2[ . 
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FIGURES 

FIG. 1. The first excitation energy Ei{a, b) — Eo{a, h) multiplied by chain length L is plotted 
v.s. 1/lnL, L = 4, 6, • • • ,60, for different values of the anisotropy a = ±0.5 and impurity bond 
strength h = 0,1,0.1. 

FIG. 2. The scaled low energy levels 62(0,6) and 63(0,6) of Eq.(2) for attractive interaction 
(a = -0.5) are plotted v.s. 1/lnL, L = 4, 6, • • • , 50. e2(-0.5,0), e3(-0.5,0) is for OBC (6 = 0), 
e2(-0.5,l), 63(-0.5,l) for PBC (6= 1), and e2(-0.5, 0.1), 63(-0.5,0.1) for impurity chains. 

FIG. 3. The scaled low energy levels 62(0,6) and 63(0,6) of Eq.(2) for repulsive interaction 
(o = 0.5) is plotted v.s. 1/lnL, L = 4,6, ••• ,40. 62(0.5,0), 63(0.5,0) is for OBC (6 = 0), 62(0.5,1), 
63(0.5, 1) for PBC (6 = 1), and 62(0.5,0.1), 63(0.5,0.1) for impurity chains. 

FIG. 4. The exponents a(o, 6) and /3(a, 6) of Eq.(lO) are plotted for both repulsive and 
attractive interactions o = ±0.5 and impurity bond strength 6 = 0.1,0.9 as a function of 
-2/[lnL + ln(L + 2)] with L = 4, 6, • • • , 48. The fitting lines indicate that a(-0.5,0.1) and the 
/3(0.5,0.1) flow to zero, while q;(0.5,0.1) and /3(-0.5,0.1) flow to 1/16. 
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Fig.2 Qin, Michele, and Yu 
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